Formation and Transformation of Vector-Solitons in Two-Species Bose-Einstein 

Condensates with Tunable Interaction 
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Under a unified theory we investigate the formation of various types of vector-solitons in two- 
species Bose-Einstein condensates with arbitrary scattering lengths. We then show that by tuning 
the interaction parameter via Feshbach resonance, transformation between diff'erent types of vector 
solitons is possible. Our results open up new ways in the quantum control of multi-species Bose- 
Einstein condensates. 
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Since the first realization of Bosc-Einstcin condensate 
(BEG), a tremendous amount of research has taken place 
in this interdisciplinary field P, Q ■ One specific topic of 
wide interest is multi-component BECs, which possess 
complicated quantum phases [sl and properties not seen 
in individual components. Theoretical [j, [B| and ex- 
perimental [3, H, Q studies have shown that inter-species 
interaction plays crucial roles in these systems. Recently, 
two-species BECs with tunable interaction [l^, [HI were 
successfully produced . This important progress moti- 
vates further explorations of two-species BECs. 

Solitons represent a highly nonlinear wave phe- 
nomenon with unique propagation features and have at- 
tracted great interests from many different fields, e.g., 
nonlinear optics [T^]- Because the nonlinear Schrodinger 
equation governing matter-wave solitons 0, [l^ [3, EBl 
is similar to those for solitons in other contexts, well- 
established understandings of solitons can be applied to 
atomic BECs to achieve better control of matter waves. 
However, BEC systems possess many unique features of 
their own and make it possible to study soliton phenom- 
ena under circumstances that cannot be easily realized 
in other fields. 

In this paper, we provide a unified theory to investigate 
various types of vector solitons in two-species conden- 
sates with arbitrary scattering lengths. Remarkably, we 
find that by tuning the inter-species interaction via the 
well-established Feshbach resonance technique [l^, dif- 
ferent types of vector-solitons can be realized and even 
transformed to each other. In addition, these vector 
solitons can be dynamically stable in the presence of 
soft trapping potentials. Our theoretical work opens up 
new possibilities in the quantum control of multi-species 
BECs and other related complex systems. 

The mean-field dynamics of a two-species BEC is gov- 



erned by the following equations 0, 0, H, 0] : 
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where the condensate wave functions are normalized by 
particle numbers Ni = J \ipi\'^d^r, Uu = Airfi^aii/mi and 
Ui2 — U21 = 2iTh'^ai2/m represent intra- and inter- 
species interaction strengths respectively, with atj be- 
ing the corresponding scattering lengths and m the re- 
duced mass. The trapping potentials are assumed to 
be Vi — mi[uLif^x^ + ^i^iv"^ + ^^)]/'^- Further assuming 
uji± ^ ujix such that the transverse motion of the con- 
densates are frozen to the ground state of the transverse 
harmonic trapping potential, the system becomes quasi- 
one-dimensional. Integrating out the transverse coordi- 
nates, the resulting equations for the axial wave functions 
"01,2(2^) in dimensionless form can be written as 
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6l2|V'2|1V'l, (1) 



^ = [-^^ + + b2,m' + b22\U]i^2, (2) 



where we have chosen ^Jh/ [miUJw^ and 27: /ujn_ to be the 
units for length and time, respectively; and ■01^2 is nor- 
malized such that / \%j}i\'^dx = 1 and / \'ip2\'^dx — N2/N1. 
Other parameters in Eqs. ((T)) and ^ are defined 
as: bii = 2aiiNi, 612 = 621 = 2miai2iVi/[(l + 7)771], 
622 = 2a22miNij/m2, 7 = uj2±/'^i±' ^1 = ^ix/(^i±, 
A2 = W2a;/wij_, and K = TOi/m2. 

A unified treatment of different vector-soliton solu- 
tions - When the longitudinal trapping potential is ne- 
glected {\i = 0), Eqs. Il]) and ([2]) are perfectly inte- 
grable under the condition k = 1 (i.e., mi = TO2) and 
bii = bi2 = 621 — 622 [3, allowing for a general pro- 
cedure to construct two-component vector-soliton solu- 
tions in the form of "dark-dark" \it\, "bright-dark" [11] 
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and "bright-bright" sohtons [l9|. When integrabihty is 
destroyed, previous studies show that distorted versions 
of soliton solutions exist, but closed form can only be 
given for special cases and often for one particular type 
of vector soliton 20] . Here we show that, even when the 
above-mentioned integrabihty condition is violated (e.g., 
for arbitrary interaction strengths bij and mass ratio k), 
there still exists in general a specific class of exact solu- 
tions that include all types of vector solitons, so long as 
bi2 ^11^22- The type of vector soliton is determined by 
the values of bij and k. 

The vector-soliton solution can be found by inserting 
an appropriate ansatz (see below) into Eqs. ^ and ([2]). 
Defining the following two quantities: 



Ci 

C2 



(622 - k6i2)/(&?2 - &II&22) 

(612 - k6ii)/(6i2 ~ &11&22) 



the conditions under which various vector-soliton solu- 
tions exist are found to be: 



Ci > 0, C2 < 
Ci > 0, 6-2 > 
Ci < 0, C2 > 
Ci < 0, C2 < 



bright - bright (BB), 
bright - dark (BD), 
dark - dark (DD), 
dark - bright (DB). 



Explicit expressions for these vector-solitons are then 
found from the ansatz we use. In particular, the BB so- 
lution (i.e., bright soliton for species 1 and bright soliton 
for species 2; analogous convention applies to all other 
solutions) is given by 
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77Vt^ sech(?72; - r]vt) e4"^+('''-^')*/2] , 



The BD solution is given by 

'4>iB = r]\/cIsech{T]x - T]vt) e'[^=^+/i*l , (3a) 
^2D = [ivy/C^/K + f]y/C^ ta,iih{rix-r]vt)y^^\ (3b) 



with /i = (77^ 



/2 - 6i2C2(77' + v^K.^) and ^ 



-622C2(?7 -I- w /k ). The DD solution is found to be 



•4>iD = [iv\/-Ci + T]\/ -C\ tanh(77a; - Tywt)] e*-^^* , 
■021) = [ivy/Ch./ n + rjy/Ch. tanh(r/a; — rjvt)] e*-^^* , 

with/i = -if -v'^(-biiCi+bi2C2l ^J^) and /2 = -fJ-V^- 
v^(— 621C1 + ^22C'2/m^)- Finally, the DB solution can be 
obtained from the BD solution ([3]) by exchanging indices 
1 and 2, and let Ci_2 — > — Ci^2- In all these cases, the 
parameter rj determines the width of the soliton and can 
be found by the normalization condition for V'1,2. The 
parameter v gives the velocity of the soliton. We note 
that the BB solution under the condition 611,622 > 
and 612 < 0, but not the other cases, has been studied by 
Perez-Garci'a and Beitia [2l|. The exact expressions for 
vector solitons found here bear an apparent resemblance 
to those for integrable systems [l6j . 
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FIG. 1: Phase diagram of vector solitons vs. inter-species 
interaction 612 for the system with 611 < and 622 > 0. As 
described in the text, the unit of bij is \/h/ (mitJix). In the 
examples shown later in the paper, the two species experience 
the same trapping frequency: u}i± = u>2± = 2n x 710Hz. The 
number of atoms in the first species (^Li) is assumed to be 
5000, which results in the dimensionless intra-species interac- 
tion strength 611 = —1.47 and 622 = 5.89. The two critical 
values of the inter-species interaction strength are therefore 
Kbit = -0.45 and b22/« = 19.35. 



Dynamics and dynamical stability - To be specific, we 
take ^Li as species 1 and ■^'^Na as species 2, with fixed 
and realistic intra-species interaction strength 611 < 0, 
622 > and variable inter-species interaction strength 
612. Figure [1] depicts the "phase diagram" for different 
regimes of vector solitons as 612 is varied. As is indi- 
cated in Fig. m the system supports BB (612 < k6ii), 
BD (k6ii < 612 < 622/K) and DD (612 > 622/K) vector 
soliton, while the condition for DB soliton can never be 
satisfied for this particular ''Li-^'^Na system. 

To study the dynamical stability of the vector-soliton 
solutions given above, we first calculate the collec- 
tive excitation frequencies by solving the corresponding 
Bogoliubov-de Gennes equation Our calculations 

show that, the excitation frequencies of the BB and BD 
vector-solitons are all real, while those of the DD vector- 
soliton contain complex values. We hence conclude that, 
for the ^Li-^'^Na system, the BB and BD vector-solitons 
are dynamically stable, while the DD vector-soliton is 
dynamically unstable. 

This stability analysis is confirmed by direct numerical 
simulations of Eqs. (H]) and Computational examples 
are shown in Fig. O In our simulation, a harmonic trap- 
ping potential is added along the longitudinal direction to 
represent a more realistic situation. The trapping poten- 
tial is weak such that its variation across the soliton scale 
is negligible. Such a soft trapping potential spatially con- 
fines the solitons without affecting their essential prop- 
erties. The upper panel of Fig. [2] shows the evolution 
of a BD vector soliton. For the initial condition, we use 
the exact solution of Eq. ([5a|) for the bright component, 
while in order to confine the dark soliton inside the trap, 
we multiply Eq. (|3bp by a Thomas-Fermi profile with an 
inverted parabolic shape to simulate the dark component 
which is a quite standard practice [T8| . The presence of 
the weak longitudinal trap causes the vector soliton to 
oscillate while maintaining the overall shape as shown in 
the upper panel of Fig. [2l By contrast, the lower panel of 
Fig. m shows that the DD soliton is dynamically unsta- 
ble — an initial DD vector-soliton is quickly destroyed 
and the two species tend to phase separate into different 
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FIG. 2: (Color online) Dynamics of vector solitons inside a 
weak longitudinal trap with Ai = A2 = 0.1. Other parameters 
are specified in Fig. [T] Units for x and time are \/h/ (miuix) 
and 27r/tJix. The left and right plots represent the longitu- 
dinal density profiles of the first and second species, respec- 
tively. Upper panel: evolution of a BD vector soliton with 
612 = 2 initially located at s = and velocity v = 0.1. Lower 
panel: evolution of a DD vector soliton with 612 = 30. For 
the system under study, the BD (DD) vector soliton is dy- 
namically stable (unstable). 



spatial regions due to the large inter-species repulsion. 
We remark that there exist alternative ways to study the 
soliton dynamics, for example, the variational approach 
used in Ref. [2^ . 

Soliton transformation via Feshbach resonance - Con- 
trolling the interaction strength via Feshbach resonace is 
a unique possibility afforded by atomic systems. Fig. [3] 
shows an example of dynamical conversion between dif- 
ferent types of vector solitons via an inter-species Fesh- 
bach resonance. Initially, we prepare a zero-velocity BD 
vector soliton by choosing a proper &12 whose value is 
tuned in later time. In the upper panel of Fig. [31 612 
is changed to a value less than nbn. For this new value 
of &12, the system supports BB, but not BD, vector soli- 
tons (see Fig. [T|) . The upper panel of Fig. [3] shows that 
the initial BD vector soliton is indeed destroyed after &12 
is changed to the new value. The system first evolves 
into two pairs of BB vector-solitons with finite relative 
velocity. These two pairs go through several collisions 
and eventually merge into one single pair. During the 
evolution, significant radiation from the second species 
is observed. According to the exact BB solution given 
above, the number of atoms in each species satisfy 

Ni _ Ci _ 622 - ^^12 
N2 -C2 ubii - bi2 ' 

Using realistic values of 611 and 622 (as detailed in the 
caption of Fig. 1), we have N1/N2 — 0.935. Numeri- 
cally, we found this ratio to be 0.92 0.95 which is in 




FIG. 3: (Color online) Conversion between different types of 
vector solitons via Feshbach resonance. The left and right 
plots represent the longitudinal density profiles of the first 
and second species, respectively. Units for x and time are 
^^/(miajix) and 2tt/uji±. Initially, a BD vector soliton is 
prepared with &12 = 2. Other parameters are the same as in 
Fig. [2] Upper panel: From t = 3 to 4, 612 is linearly varied to 
a final value of —10 and stays at that value afterwards. Lower 
panel: From t = 3 to 23, bi2 is linearly varied to a final value 
of 30 and stays at that value afterwards. 



good agreement with the theoretical value. It is sim- 
ple to check and intuitive to understand that when the 
initial BD vector-soliton breaks into two pairs of BB soli- 
tons, the relative phase of the two bright solitons of the 
first species is zero, while that of the second species is tt. 
Just like in the scalar condensate case, this leads to at- 
tractive (repulsive) soliton-soliton interaction within the 
first (second) species. We also remark that the relative 
velocity of the BB vector-soliton pairs thus generated is 
sensitive to how fast 612 is varied. In general, the faster 
612 is varied, the larger the relative velocity. If we vary 
612 much slower, the pairs (with smaller relative velocity) 
will merge together later and experience more collisions. 
This observation suggests that varying the sweeping rate 
of bi2 should help understand more about the collision 
dynamics of multiple BB vector-soliton pairs [21], . 

In the example shown in the lower panel of Fig. [31 612 
is changed to a value larger than 622 /'v, i-C, into the DD 
vector-soliton regime. Here the BD vector soliton does 
not evolve into a DD vector soliton, which is consistent 
with the fact that the DD solution is unstable. Instead, 
due to the strong inter-species repulsion, the core region 
of the initial dark soliton of species 2 is widened, and 
is no longer accurately described by a tanh function. In 
the mean time, the bright soliton of species 1 is narrowed 
and remain captivated inside the core of the dark com- 
ponent. This property may help to produce and stabilize 
bright solitons. During the evolution, we can also see that 
additional dark solitons can be generated from the core 



4 



region of species 2 (see the lower right figure in Fig. [3]). 
Of particular interest are those new dark solitons gener- 
ated at early times: They move away from the core, then 
get bounced back by the trapping potential, and finally 
get reflected by the central BD vector soliton. Recently 
similar behavior — the reflection of dark solitoi i by a BD 
vector soliton — was observed experimentally [l5| . 
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FIG. 4: (Color online) Conversion between BB and BD vector 
solitons. The left and right plots of the upper panel represent 
the longitudinal density profiles of the first and second species, 
respectively. The lower panel shows the variation of 612. Unit 
for X and 612 is y^h/{miu!i±), and unit for time is 2-k/lji±. 
Initially, a BD vector soliton is prepared with &12 = 2. Other 
parameters are the same as in Fig. (2] 

Figure 2] illustrates another example of the conversion 
between different types of vector-solitons via inter-species 



Feshbach resonance. The initial condition is the same 
as in Fig. [3] where a zero- velocity BD vector-soliton is 
prepared. The value of bi2 is then decreased and the 
system evolves to two BB vector-soliton pairs as in Fig. [31 
Then, we restore the initial value of 612 before the two 
BB vcctor-solitons merge into one. As seen in Fig.lH such 
a controlled tuning of 612 can also recover the initial BD 
vector-soliton. This represents a remarkable example of 
the coherent control over the BEG dynamics, afforded by 
the tunability of interaction strength. 

In conclusion, we put forward a unified theory to treat 
all types of vector solitons on equal footing in two-species 
BECs with arbitrary scattering length. More impor- 
tantly, we show that different types of vector solitons 
can be transformed to each other by varying the interac- 
tion strength. Besides its fundamental interest in quan- 
tum dynamics, our study may also find applications in 
the quantum control of multi-species BECs. For exam- 
ple, treating one atomic species as the target, its soliton 
properties can be controlled by the presence of a second 
species via the tuning of the interaction strength between 
them. We hope that these results can stimulate investiga- 
tion of novel phenomena in nonlinear complex systems in 
general and further studies of BEC solitons in particular. 
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